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I. INTRODUCTION 

Experimental studies on charmonia decay properties are essential to test perturbative QCD models and QCD 
based calculations. The importance of the Color Octet Mechanism (COM) for Xc.j decays has been pointed 
out for many years 0|, and theoretical predictions of two-body exclusive decays have been made based on it. 
Recently, new experimental results on Xcj exclusive decays have been reported 0, 1^ . COM predictions for 
many Xcj decays into meson pairs are in agreement with experimental values, while predictions for some decays 
into baryon pairs, for example, the branching fractions of Xcj AA, disagree with measured values. For further 
testing of the COM in the decays of the P-wave charmonia, measurements of other baryon pair decays of Xcj, 
such as Xcj ^ and are desired. 

The measurement of Xco is helpful for understanding the Helicity Selection Rule (HSR) which 

prohibits XcO decays into baryon antibaryon {BB) pairs. However, the measured branching ratios for XcO decays 
into pp and AA do not vanish, demonstrating a strong violation of HSR in charmonium decays. Measurements 
of XcO decays into other baryon anti-baryon pairs would provide additional tests of the HSR. 

In this paper, the analysis of ipi'^S) — > 72(7r+7r^)pp decays using 14 x 10^ tpi'^S) events collected at BE- 
SII/BEPC is reported. We observe Xcj S^S+ and AAtt+tt^ decays and search for Xcj KgKgpp; results 
are compared with theoretical predictions. 

II. BES DETECTOR 

This analysis is based on a sample of 14 x 10^ '0(25') events taken with the BESH detector at the BEPC 
storage ring at a center-of-mass energy of M^(^2S) with an integrated luminosity of 19.72±0.86 pb~^ jH . BES is 
a conventional solenoidal magnet detector that is described in detail in Ref. Q ; BESH is the upgraded version 
0- A 12-layer vertex chamber (VC) surrounding the beam pipe provides trigger and trajectory information. 
A forty-layer main drift chamber (MDC), located radially outside the VC, provides trajectory and energy 
loss (dE/dx) information for charged tracks over 85% of the total solid angle. The momentum resolution 
is (Jp/p — 0.017^/1 -I- p^ (p in GeV/c), and the dE/dx resolution for hadron tracks is ~ 8%. An array of 
48 scintillation counters surrounding the MDC measures the time-of-flight (TOE) of charged tracks with a 
resolution of ^ 200 ps for hadrons. Radially outside the TOE system is a 12 radiation length, lead-gas barrel 
shower counter (BSC). This measures the energies of electrons and photons over ~ 80% of the solid angle with 
an energy resolution oi ge/E — 22% /^/E (E in GeV). Outside of the solenoidal coil, which provides a 0.4 Tesla 
magnetic field over the tracking volume, is the iron flux return which is instrumented with three double layers 
of counters that identify muons with momentum greater than 0.5 GeV/c. 

III. MONTE-CARLO SIMULATION 

A GEANT3 based Monte Carlo simulation program, SIMBES which simulates the detector response, 
including interactions of secondary particles in the detector materials, is used to determine detection efficiencies 
and mass resolutions, as well as to optimize selection criteria and estimate backgrounds. Under the assumption 
of a pure El transition, the distribution of polar angle 9 of the outgoing photon in ■0(25') — > 7Xcj decays is 
given by 1 -I- acos^6' with a = 1, —1/3, and 1/13 for J = 0, 1, and 2, respectively @. Angular distributions 
of daughter particles for the sequential decays Xc,, S Att are simulated based on the transitional 

amplitude information method Angular distributions of daughter particles from the other multi-body 

decays are generated isotropically in the center-of-mass of ip{2S) or Xcj- 

IV. EVENT SELECTION 

The selection criteria described below are similar to those used in previous BES analyses 0. 

A. Photon identification 

A shower cluster in the BSC is considered to be a photon candidate, if it has an energy deposit of more than 
50 MeV, the angle between the nearest charged track and the cluster is greater than 12°, the first hit is in the 
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beginning six radiation lengths, and the difference between the angle of the cluster development direction in 
the BSC and the photon emission direction is less than 37°. 



B. Charged particle identification 



Each charged track is required to have a good helix fit and a polar angle 9 that satisfies \cos9\ < 0.8. The 
TOF and dE/dx measurements of the charged track are used to calculate Xp/^) values and the corresponding 
confidence levels for the hypotheses that the particle is a pion, kaon, or proton. 



C. Event selection criteria 



Candidate events are required to satisfy the following selection criteria: 

1. The number of charged tracks is required to be six with net charge zero, and the number of photon 
candidates must be less than four. 

2. The proton and antiproton must be identified using particle identification; their particle identification 
confidence levels must be greater than 1%. 

3. Four-constraint (4-C) kinematic fits to the i^{2S) —> j2(Tr^TT^)pp hypothesis are performed using each 
photon candidate. The combination with the minimum combined Xcomb ~ S"'^ iXdE/dx~^^TOF)i~^X4c 
is selected, and the 4-C fit confidence level for the selected combination must be greater than 1%. 

4. To remove the background channels: i^{2S) 72(7r+7r~)i4r+i4r", 73(7r+7r~), and 2(7r+7r~)pp, Xcomb 
signal channel is required to be less than those of the background channels, i.e., x^g^i,{signal) < Xcombi^9)- 



V. DATA ANALYSIS 



Xcj 



The sequential decays S A7r,A pn are used to reconstruct Xcj The pTT~ invariant mass 

distribution for candidate events is shown in Figure Qfa), where a clear A signal is seen. A A signal appears 
similarly in the pTT~^ invariant mass distribution. The A mass resolution, obtained from MC simulation, is 
5.5 MeV/c^. After requiring that A and A satisfy |mA — mp7r| < 11 McV/c^ (Sct), the scatter plot of toa^- 
versus mj^^+ is given in Figure ^^b). In Figure ^c), the ttt-att- distribution is shown after the requirement 
l'^A7r+ ~ '^hI < 17.5 MeV/c^ (2c), where the S~ appears clearly. After requiring that Avr^ and A7r+ satisfy 
\Matt — M-e\ < 17.5 MeV/c^, the TO2-h+ distribution is plotted in Figure^d), and the normalized background 
distribution determined from the A A sideband region is also shown in the plot. 

To determine the detection efficiencies, the angular distributions of the decay particles for the sequential 
decays ^^{23) —> IXcjtXcj ~^ S~S"'",5 — > A7r,A — > pi: are simulated based on the transitional amplitude 
information (TAI) method 0^ The detection efficiencies are (2.3 ± 0.1)%, (2.6 ± 0.1)%, and (2.3 ± 0.1)% for 
ip{2S) sequential decays to 3~S+ via Xco, Xci, and Xc2, respectively. 

The main backgrounds from channels with S or A production, including ^1^{2S) E!~S+, ^'(25) "fXcj 
7l](1385)fi(1385),V'(25') ^ 7Xc,, ^ 7S"fio, and V(25) -fXc., -> 1^°^° are determined by Monte-Carlo 
simulation. By using the branching ratio of ip(2S) — > S~S+ measured in |llj | and naively assuming Xcj ~* 
E(1385)i;(1385),E''E°,S°S'' have the same branching ratio as Xc., ^ AA, one obtains 0.03, 0.01, and 0.02 
background events for XcO, Xcii and Xc2, respectively. 

Figure 121 shows the fit to the m^-2+ distribution using a signal shape obtained from MC simulation. One 
obtains 6.4 ± 3.2 events for XcO ~^ S~S+. Within twice the mass resolution of the Xci and Xc27 one event 
is respectively observed. Upper limits at the 90% Confidence Level (CL) for Xci and Xc2 S~S+ 

are evaluated with POLE 12] including systematic uncertainties, and the results are listed in Table U The 
branching ratios and upper limits (A^''*^ —> Nipper) are estimated with: 

N''"''/e 

B[xcJ - ^ = ] - ^^^^^^5[^(25) _ ^^,j]B[^ -> ^A]2S[A ^ pirY ' 
where e denotes the detection efficiency and TV^s = N°'^^ — N^^ denotes the number of signal events observed. 
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Fig. 1: (a) The pn invariant mass distribution, (b) Scatter plot of m^^- versus m^.^+, where 
the small square corresponds to the H~E"^ signal region and the larger closed area to the H~H"'' 
sideband region, (c) Att^ mass distribution for |m^^+ — ms] < 17.5 MeV/c^. (d) The 771=-=+ 
invariant mass distribution, where the shaded histogram corresponds to the contribution from 
the H~H^ sideband region (normalized). 



TABLE I: Summary of numbers used in the branching ratio calculation, where the branching ratio for ^{2S) 7Xc./ 
are taken from Ref. jlSj 



quantity 




XcO 




Xc2 






6.4 ±3.2 


1.0 


1.0 


^ upper 




12.4 


4.6 


4.6 


e{%) 




2.3 ±0.1 


2.6 ±0.1 


2.3 ±0.1 


N^i2S) (xlO«) 




14.0 ± 0.6 


14.0 ± 0.6 


14.0 ± 0.6 


B{i,i2S) ^ 7Xc,7 


) (%) 


9.2 ±0.5 


9.1 ±0.6 


9.3 ±0.6 


B{E ^ Att) (%) 




99.9 ±0.1 


99.9 ±0.1 


99.9 ±0.1 


B{A ^ pn) (%) 




63.9 ±0.5 


63.9 ±0.5 


63.9 ±0.5 


BixcJ - 


X 10"" 


5.3 ±2.7 ±0.9 


< 3.4 (90% C.L.) 


< 3.7 (90% C.L.) 



< 10.3 (90% C.L.) 
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Fig. 2: Fit to the mass spectrum for XcO E~E'^ with the signal shape obtained from 

Monte-Carlo simulation. The histogram corresponds to the data, and the shaded histogram to 
the normalized backgrounds obtained from Monte-Carlo simulation. 



B. Xcj AAtt+tt" 

The A and A candidates are reconstructed the same as in the above analysis. The Xcj candidates 
are excluded. Figure|3ta) shows the scatter plot of nip^- versus ?Tip^+ for candidate events. The dense cluster 
in the square, with side twice the A mass resolution, corresponds to the AA signal. The AA sideband region 
is shown in the same plot. Figure I3b) shows the distribution of ni/^f^^^ (full histogram) and A A sideband 
background (shaded histogram). 

To determine the detection efficiency, sequential decays via Xc,, AAtt+tt", A — > pn are simulated assuming 
a pure phase space distribution. After applying the same selection criteria, one obtains detection efficiencies of 
(0.93 ± 0.04)%, (1.14 ± 0.04)%, and (1.11 ± 0.04)% for XcO, Xci, and Xc2 decays, jespectively. The decrease in 
detection efficiency, compared with the two-body sequential decay via Xcj ~^ S~S^", is due to the difference in 
the momentum distributions of the decay particles in Xcj AAtt+tt". 

Potential background contaminations, including 'ip{2S) — > 'fXcj 7p*'A++A^^ 72(7r+7r^)pp, ip{2S) 
7r"*'7r~ J/'0, then J/ip —> jijc — > jpn^A ,^p^pp, and J/%1: — > rjpp —> jTr^w^pp, are studied by Monte-Carlo 
simulation, and it is found that these background contaminations are effectively removed by the AA sideband 
background subtraction. 

Figure 0] shows the mj^]i^^+^- distribution after subtracting the AA sideband background. Breit-Wigner 
functions, convoluted with Gaussian mass resolutions, for the XcOiXcOj and Xc2 resonances and a polynomial 
background shape are used to fit the mAA^+^- distribution from 3.3 to 3.65 GeV/c^. The number of signal 
events obtained are 71**'^ = 9.6 ± 5.4, 0.2 ± 3.2, and 10.4 ± 5.7 for XcO, Xa, and Xc2, respectively. The upper 
limits at the 90% C.L. are estimated with POLE including systematic uncertainties and also listed in Table Ull 



TABLE II: Summary of numbers used in the branching ratio calculation for Xcj ^ AA7r''"7r 



quantity 




Xcl 


Xc2 


n ^ 


9.6 ±5.4 




10.4 ±5.7 


Supper 


19.6 


8.8 


20.4 


6 (xio-^) 


9.3 ±0.4 


11.4 ± 0.4 


11.1 ±0.4 


N^(2S) (xlO«) 


14 ±0.6 


14 ±0.6 


14 ±0.6 


B{^{2S) ^ 7Xcj) (%) 


9.2 ±0.5 


9.1 ±0.6 


9.3 ±0.6 


B(A ^ p-k) (%) 


63.9 ±0.5 


63.9 ±0.5 


63.9 ±0.5 


b{xc.J TT + TT-AA) X 10"'' 


2.0 ± 1.1 ± 0.4 (2.5a) 




1.8 ± 1.0 ±0.3 (2.5cr) 


upper-limit xlO"^ (90% C.L.) 


< 4.0 


< 1.5 


< 3.5 
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Fig. 3: (a) Scatter plot of m^^- versus mp^+ for candidate events, where the square box cor- 
responds to the AA signal region, and the large box beside it corresponds to the AA sideband 
region, (b) The my^^TTTr invariant mass distribution of events in the AA region (full histogram) 
and normalized AA sideband background (shaded histogram). 
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Fig. 4: Fit to the mAA7r+7r- distribution after subtracting AA sideband background. The squares 
with error bars are data, the dashed curve is a second order polynomial to represent background, 
and the solid curve is the fit. 



The branching ratios and upper limits [N''^^^ Nipper) are estimated using: 



where n''*^ denotes the number of signal events. The branching ratios and upper limits are listed in Table Hll 
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C. xcj ^ K%K%pp 

Kg candidates are selected by choosing the tt^tt^ combinations with the minimum value of A(ii,i2, ji, j2) = 
- rriiio)'^ + {m^+(^2)7T-U2) - "^ifg)^ and requiring the ifg candidates to satisfy |M^o - M^+^- \ < 
15 MeV/c^ {2a). The candidates for Xc, ^ AA7r+ TT are excluded. 

The scatter plot of m^+^- versus m^+^~ is shown in FigureEta), where the events in the central square are 
KgKg signal candidates and the four surrounding squares are the KgKg sideband regions. Figure l^b) shows 
the mTTTT invariant mass distribution. In the Kg mass region, no significant signal is observed. After applying 
the selection criteria for KgKg, the distribution of mj^oj^Opp is given in Figure [Sfc). Within selection regions 
which are twice the Xcj mass resolution, the observed numbers of events are 2, 0, and 2 for Xco, Xcii E^nd Xc2, 
respectively. FigureEfd) shows the background from the KgKg sideband background region. The Kg sideband 
regions are defined as 445.1 MeV/c^ < m^+^- < 475.1 MeV/c^ and 520.1 MeV/c^ < m^+^- < 550.1 MeV/c^, 
and the mass of other two pions is required to satisfy \Mj^o— M^+^- \ < 15 MeV/c^. The numbers of background 
events observed are 1.3, 1.5, and 0.5 for XcO, Xa, and Xc2 respectively. The upper limits at the 90% C.L. are 
estimated with POLE including systematic uncertainties, and the results are listed in Table ITTTl 




> 

0) 



^ 20 



.25 0.525 0.8 

m(7T*(1)7T"(2)) (GeV/c') 




0.2 0.5 0.. 

m{-n*n') (GeV/c') 



1 



> 

°0.75 



: 0.5 
j 

0.25 





3.4 3.6 
m(KsKspp) (GeV/c") 




3.4 

m(KsKspp) (GeV/c') 



3.6 



Fig. 5: (a) The scatter plot of m^+^^ versus m^+^- . (b) The m^+^- invariant mass distribution, 
(c) The rnifOf^Opp invariant mass distribution, (d) The JTi/fgifgpp invariant mass distribution 
(normalized) for events in the Kg Kg sideband regions. 



The detection efficiencies are (0.87± 0.06)%, (0.92 ± 0.06)%, and (0.90± 0.06)% for XcO, Xa, and Xc2 decays, 
respectively. Upper limits are estimated using 

B[x.j K^KIpp] = (3) 

N^{2S)B[ip(2S) -fXcj]B[K^ vr+TT ]^ 

where e denotes the detection efficiency. The upper limits for Xcj KgKgpp are listed in Table ITm 
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TABLE III: Summary of numbers used in the branching ratios calculation of Xcj KgKgpp. 



quantity XcO Xci Xc2 



''Upper 



4.7 2.5 4.4 



e(xlO-') 8.7±0.6 9.2 ± 0.6 9.0 ± 0.6 

N^(2S) (xlO^) 14.0 ±0.6 14.0 ±0.6 14.0 ± 0.6 

B(V'(2S') ^ 7Xcj) (%) 9.2 ±0.5 9.1 ±0.6 9.3 ± 0.6 
bIk^ T^^T^~) (%) 69.0 ± 0.1 69.0 ± 0.1 69.0 ± 0.1 

BiXcj ^ K"sK"spp) X 10-* 

upper-limit (90% C.L.) < 8.8 < 4.5 < 7.9 



VI. SYSTEMATIC ERRORS 



Systematic errors on the branching fractions mainly originate from the MC statistics, the track error matrix, 
the kinematic fit, particle identification, the photon efficiency, the uncertainty of the branching fractions of the in- 
termediate states (from PDG or published papers), the total number of ip{2S) events, the fitting, the uncertainty 
of the angular distribution of the S in Xcj decays, and efficiency corrections from Xc, I]-(1385)A7r+ ± c.c. 

1. As studied by a previous BES analysis [l| using clean channels like J/4> AA and ■0' ^ tt+tt J/ip ^ 
7r+7r~/i+/i-, it is found that the tracking efficiency for MC simulation agrees with data within (1-2)% for 
each charged track. Hence, we take 12% as the systematic error for events with six charged tracks. 

2. The photon detection efficiency was studied using different methods with J/tp — + t:^tt^tt^ events [T5| . 
and the difference between data and MC simulation is about 2% for each photon. We take 2% as the 
systematic error for decays with one photon. 

3. A kinematic fit is applied to the channels analyzed, and a probability > 0.01 is required. The consistency 
between data and Monte Carlo efficiencies for the kinematic fit depends on whether the track fitting error 
matrices are consistent for data and MC simulation. From earlier studies, we take 4% as the systematic 
error for kinematic fitting |l6l |. 

4. The p and p candidates are identified with the requirement that the particle identification confidence level 
be greater than 0.01. Comparing data and MC data samples for the decay tt+tt" J/t/i, J/i/j 
TT^pp, the efficiency difference is about 2.8% for one particle identified as a proton or antiproton. Here 
5.6% is taken as the systematic error for pp identification. 

5. To estimate the background uncertainty for Xcj ~^ S~S"^, background MC samples and S^S"*" sideband 
backgrounds are used. The maximum difference between them is about 5%, which is taken as the system- 
atic error. For Xcj AAtt+tt^ , the background uncertainty includes the uncertainty in the background 
shape and fitting interval; the differences found using different shapes and intervals are 5.25%, 3.57%, and 
2.85% for Xco, Xcii and Xc2 decays, respectively. For Xcj K^Kgpp decays, the maximum difference in 
estimation of the upper limits between different choices of KgKg sideband regions is about 4%, which is 
taken as the systematic error. 

6. To determine the detection efficiency for Xcj the TAI method is used to account for the 
angular distribution of the 2. The uncertainty of the parameters used introduce 4% and 5% differences 
in efficiencies of the Xci and Xc2 channels, respectively, which are taken as the systematic error for the 
uncertainty of the S angular distribution. 

7. In the analysis of Xcj AAtt+tt^, the S~S+ intermediate state is excluded, and no significant signals 
for decays Xc,, ~> S+(1385)S~(1385), A(1405)A(1405) are observed. The intermediate state corrections to 
the detection efficiency are evaluated by considering only the observable decay Xcj E~(1385)A7r+±c.c.. 
The branching fractions are roughly estimated using data. The efficiency differences are evaluated by: 

Se _ |ei - 62! B{xcj ^ S-(1385)A7r+ + c.c.)B{^- Avr") 

ej" ~ 7i B{xcj S-(1385)A7r+ + c.c.)B(E- ^ Att-) + B{xcj ^ AAtt+tt-) ' 

where ei and £2 are the detection efficiencies for ^p{2S) radiative decays via Xci ~* AAtt+tt^ and Xcj ~* 
E"(1385)A7r+ ± c.c, i;(1385) -> Att, respectively. We obtain the efficiency difference 6e/ei = 7.1%, 6.0%, 
and 2.6% for decays of XcOj Xci, and Xc2, respectively, which is treated as one of the systematic errors. 
Total systematic errors are given in Table IIVI 
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TABLE IV: Summary of systematic errors (%). 
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XcO 


Xc2 


iViJJU tracK 


ion 


ion 


ion 


ion 
iz.u 


12.0 


12.0 


ion 
iz.u 


ion 
iz.u 


ion 
iz.u 


Photon efficiency 




z.U 


z.U 


n 
z.U 


2.0 


2.0 


o n 
z.U 


o n 
z.U 


n 
Z.U 


4C-fit 






4 


4 n 


4.0 


4.0 


4 n 


4 n 


4 n 


Number of ^^(25) 


4.0 


4.0 


4.0 


4.0 


4.0 


4.0 


4.0 


4.0 


4.0 


particle ID 


5.6 


5.6 


5.6 


5.6 


5.6 


5.6 


5.6 


5.6 


5.6 


B{K pTv) 


0.8 


0.8 


0.8 


0.8 


0.8 


0.8 


0.8 


0.8 


0.8 


B{^{2S) jXcj) 


5.1 


6.4 


6.7 


5.1 


6.4 


6.7 


5.1 


6.4 


6.7 


Background uncertainty 


5.0 


5.0 


5.0 


5.3 


3.6 


2.9 


4.0 


4.0 


4.0 


MC statistics 


2.0 


1.9 


2.1 


3.5 


3.0 


3.0 


5.0 


5.0 


5.0 


Angular distribution of 5~E"'" 




4.0 


5.0 














Intermediate state S(1385) 








7.1 


6.0 


2.6 








Total 


16.3 


17.2 


17.6 


18.1 


17.6 


16.7 


16.7 


17.1 


17.3 



VII. RESULTS AND DISCUSSION 



Using the 14 x 10^ V'(2S') events accumulated at BESII, the analysis of Xcj E^E^ , AAtt^tt^ , KgKgpp 
is carried out. The measured branching fractions or upper limits are summarized in Table along with 
some theoretical predictions. Theoretically, the quark creation model (QCM) predicts B{xcO — 
(2.3 ± 0.7) X 10^^, which is consistent with the experimental value within la. For Xci and Xc2 decays into 
H~S+, the measured upper limits cover both the COM and QCM predictions. 

Within 1.8cr the branching fraction of XcO S~S+ does not vanish. For further testing of the violation of 
the HSR in this decay, higher accuracy measurements are required. 



TABLE V: The comparison of the branching fractions or upper limits for Xcj — > H^H+, AAtt+tt^ , 
and KgKgpp between experimental values and theoretical predictions. The COM predictions 
are from Ref. and the quark creation model (QCM) predictions arc from Ref. |q^ . 



Decay modes 


Branching ratios 


Theoretical predictions 






COM QCM 


XcO S"H+ 

Xcl H"H+ 
Xc2 H"H+ 


(5.3 ±2.7 ±0.9) X 10"'' 

or < 10.3 X 10"* (90% C.L.) 

< 3.4 X 10"'' (90% C.L.) 

< 3.7 X 10"* (90% C.L.) 


(2.3 ± 0.7) X 10"" 

2.4 X 10"^ 

3.4 X 10"^ (4.8 ±2.1) X lO"'^ 


XcO — » AAtt+tt" 


(2.0 ± 1.1 ±0.4) X 10"^ (2.5cr) 




Xcl AAtt+tt" 
Xc2 AAtt+tt" 


or < 4.0 X 10"^ (90% C.L.) 
< 1.5 X 10"^ (90% C.L.) 
(1.8 ± 1.0 ±0.3) X 10"^ (2.5cr) 
or < 3.5 X 10-^(90% C.L.) 




XcO ^ K'^sK"sPP 
Xcl ^ K°K°pp 
Xc2 ^ K''sK°spp 


< 8.8 X 10"* (90% C.L.) 

< 4.5 X 10"* (90% C.L.) 

< 7.9 X 10"* (90% C.L.) 
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